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The low-energy theory of interacting electrons on graphene's two-dimensional honeycomb lattice 
is derived and discussed. In particular, the Hubbard model in the large-A'' limit is shown to have 
a semi-metal - antiferromagnetic insulator quantum critical point in the universality class of the 
Gross-Neveu model. The same equivalence is conjectured to hold in the physical case N = 2, and 
its consequences for various physical quantities are examined. The effects of the long-range Coulomb 
interaction and of the magnetic field are discussed. 



A graphite monolayer, or graphene, emerged recently 
as the new frontier in physics of electronic systems 
with reduced dimensionality Such two-dimensional, 
or quasi-two-dimensional systems have led to some of 
the most startling discoveries in the condensed matter 
physics in the recent past, the quantum Hall effects and 
the metal-insulator transitions in silicon-MOSFETS and 
Ga-As heterostructures, and the high-temperature super- 
conductivity in cuprates being prime examples. What 
makes graphene qualitatively new is its semi-metallic na- 
ture with low-energy quasiparticles behaving as 'relativis- 
tic' Dirac spinors over a good portion of the conducting 
band. The spinor structure is a general consequence of 
the bipartite nature of the honeycomb lattice . Indeed, 
recently observed quantization rules for the Hall conduc- 
tivity 0] niay be understood as a direct consequence of 
the Dirac nature of its low-energy spectrum Q, [3- 

The relativistic spectrum and the concomitant lin- 
early vanishing density of states at the Fermi level, simi- 
larly as in the superconducting state of cuprates, provide 
graphene's quasiparticles with an additional protection 
against the effects of interactions. Nevertheless, a suf- 
ficiently strong repulsion is expected to turn the semi- 
metallic state into a gapped insulator, possibly breaking 
the translational and/or the rotational symmetry in the 
process. Within the simplest interacting theory defined 
by the Hubbard model there is convincing numerical ev- 
idence for the quantum phase transition at a large Hub- 
bard U into an antiferromagnet (AF) @. On the other 
hand, the long-range Coulomb interaction remains un- 
screened in the semi- metal (SM) 0, and has been ar- 
gued to favor the charge-density-wave (CDW) at strong 
coupling . The competition between different instabil- 
ities, the universality class, or even the order, of the SM 
- insulator transition, and the interplay of interactions 
with the Landau quantization in the external magnetic 
field present some of the basic open problems. Although 
graphene in its natural state may not be near a criti- 
cal point 0, one can conceive mechanical deformations 
that would pull it deeper into the strong-coupling regime 
[Tol |. Finally, the outcome of the competition between 
different interactions should have consequences for the 
selection of the ground state in the magnetic field, even 
at weak coupling [TH . 

In the present communication some of these issues are 
addressed by considering the half-filled Hubbard model 



on a honeycomb lattice, complemented with the addi- 
tional long-range Coulomb interaction between electrons. 
The analysis is based on a useful decomposition of Hub- 
bard's on-site interaction on a bipartite lattice into a 
sum of squares of average and staggered densities, and 
average and staggered magnetizations. The long-range 
part of the Coulomb interaction may be represented by a 
massless scalar gauge field, whereas its main effect on the 
lattice scale is to provide the repulsion between nearest- 
neighbors. When prepared like this, in the continuum 
limit such an extended Hubbard model on a honeycomb 
lattice maps onto a 2 + I-dimensional field theory of 
Dirac fermions, with nine different couplings. Its ap- 
parent complexity notwithstanding, when generalized to 
a large number of fermion flavors N, the theory admits a 
simple SM - AF critical point of the Gross-Neveu variety 
|l2l|. Coulomb interaction is marginally irrelevant at the 
critical point. Assuming that the equivalence with the 
Gross-Neveu model persists down to the physical case 
of = 2, I infer the values of the critical exponents in 
the original Hubbard model. A more general phase dia- 
gram, and the implications of these results for graphene 
are discussed. 

The extended Hubbard model will be defined by the 
Hamiltonian H = Hq + Hi where 



Hn = -t 



A,i,<7=±l 



ul{A)v„{A + h) + H.c. 



(1) 



Hi = 



^ 2 At:\X-Y\ 



X.Y,(y.<y' 

(2) 

The sites A denote one triangular sublattice of the hexag- 
onal lattice, generated by linear combinations of the ba- 
sis vectors ai = (-s/S, — l)(a/2), 0,2 = (0, a). The second 
sublattice is then at _B = A + 6, with the vector h be- 
ing either hi = {l/VS, l)(a/2), h = (l/VS, -l)(a/2), or 
63 = (— a/-\/3, 0). a is the lattice spacing. Neutralizing 
background is assumed, as usual. 

The doubly degenerate spectrum of Hq at E{k) = 
±i|^-exp[A: • bi]\ becomes linear and isotropic in the 
vicinity of two non-equivalent points at the edges of the 
Brillouin zone at ±K, with i? = (1, l/V3)(27r/aA/3) 0. 
Retaining only the Fourier components near ±K one can 
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write the quantum-mechanical action corresponding to 
i?o at low energies as = J^^^"^ drclxLo, with the free 
Lagrangian Lq defined as 



E 

cr=±l 



(3) 



and 



dq 



(27ra)2 

'"liK + q, UJn), uli-K + q, uj„), vl{-K + g, w„)), 

where it was convenient to rotate the reference frame so 
that qx = q ■ KjK and qy — {K x q) x K/K^, and set 
h = ks = vf = where vp = ta\/3/2 is the Fermi ve- 
locity. Choosing 70 = -^2 ® implies 71 = (8) Cy and 
72 = I2 ® cTx, with /2 as the 2x2 unit matrix, and a as 
the Pauli matrices. A « 1/a is the ultraviolet cutoff over 
which the linear approximation for the dispersion holds. 
The summation convention is adopted hereafter, but only 
over repeated space-time indices. Besides the 'relativis- 
tic' invariance, Lq also exhibits a global invariance under 
the U (4), generated by {h, ff} (Xi {/, 73, 75, 735}, where / 
is the 4x4 unit matrix, 73 = ax (E) <Ty, 75 = Uy ® Oy, 
and 735 = 17375. This is similar to the emergent 'chiral' 
symmetry of a d-wave superconductor [T^ . 

Generalizing slightly Hamman's decomposition, the 
first term in Hi can also be rewritten exactly as 

I Y^[i<^) + <A + b)f + {n{A) n{A + b)f (5) 

A 

-{m{A) + m{A + b)f - {m{A) - ■m{A + h)Y] 

where n{A),m{A) = u\{A)u+{A) ± v}_{A)u^{A), are 
the particle number and the magnetization at the sites 
A. Variables at the second sublattice are analogously 
defined in terms of v„{B). 

Defining the two slow components of the fields as 



dk 



ik-a 



l\k±K\<A (27r)2 
with r — u,v, the Dirac field becomes 



Taik^r), 



(6) 



(7) 



At low energies one may then approximate 

r„(x,T) « ri(f,T) + r2(f,T), (8) 
so that the spin densities on the two sublattices become 

rl{x,T)r„{x,T) K (9) 



-*<,(f,T)(/2 



'<Jx) ® {a^±l2)'^a{x, r). 
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FIG. 1: The large-A'^ flow diagram in tlie attractive plane 
e = gd = gj = Qx = 0- A, B, C, and G are the AF, bicritical, 
CDW, and the Gaussian fixed points, respectively. A and 
C are in the Gross-Neveu universality class. The extended 
Hubbard model in Eq. (2) defines the (dashed) line of initial 
conditions Qc = {U — V)/AUc and ga ~ —U/'lUc, with V as 
a fixed nearest-neighbor repulsion. Inset: the resulting phase 
diagram. 



with the plus sign for r = u and the minus for r = v. 

The notation is now in place to write the low-energy 
theory of the extended Hubbard model. In the continuum 
limit (a — > 0) the Lagrangian becomes 



L = L 



-^ao 



E 



E 

x—d.c^f,a 



with 



Lx = gxiY^ Wx,a'^aMx'^af 4 

(T 

gx i^Wx^a'^aMxja^.'i'a)" 
fi—3.,5 cr±l 



(10) 



(11) 



and Wd,a 
and Air = 



Mr,. = 



= 1, Wf,c = Wa^a = 0-, Af<j = Mf = 7o 

/. The short-range couplings are ga — 



-2~gd ~ e'/AK ^ {U + V)a'/8, g, ^ ~2gr ^ {U - V)ay8, 
gf ^ ga ~ ^'^gf = ^'^ga = —Ua^/8. d and c couplings 
correspond to the first (average density) and the second 
(staggered density) , whereas / and a couplings represent 
the third (magnetization) and the fourth (staggered mag- 
netization) terms in Eq. (5). The Coulomb interaction 
is represented by: 1) the intra-unit-cell, nearest-neighbor 



repulsion V 



and the 2K Fourier compo- 



nent e'^ /2K, and 2) its long-range part, which is recov- 
ered upon Gaussian integration over the scalar gauge- 
field qq. |V| should be understood as \q\ in Fourier space 
Whereas such a separation would be exact for an 
infinitely long-ranged interaction, it is only an approxi- 
mation for the Coulomb interaction. 

The usual power counting implies that all short-range 
interactions in L are irrelevant, and that the charge e 
is a marginal coupling at the non-interacting fixed point 
gx = gx = s = (i- Any critical point would therefore 
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have to lie in the strong-coupUng regime. To exert some 
control over it we may deform the Lagrangian from two 
to N flavors of the Dirac fields as follows: 

N/2 N 
o-=l a=(N/2) + l 

and ^g^/N, 2g^/N, 2e'^/N. The inte- 

gration over the Fourier components with A/b < q < A 
and —oo<uj<oo renormalizes then the short-range 
couplings at T = as 

P.^^^-9.-C^9l + 0{l/N), (13) 

- ^ = -5. + 2~gl + 0(l/iV), (14) 

with Cc.a ~ 4, Cdj = 0, and with the couplings rescaled 
as gA/ir — > 5. To the leading order in 1/N /3-functions 
for different interactions thus do not mix Since the 
model when iV = 00 is exactly solvable by the saddle- 
point method, the leading order /3-functions may also be 
understood as guaranteing that the solution is cutoff in- 
dependent . Non-analyticity of the inverse gauge- field 
propagator and the gauge invariance of L also dictate 
that 

(15) 

with z as the dynamical exponent, exactly |il4j . Rela- 
tivistic invariance of L is broken when e ^ 0, and con- 
sequently z 7^ 1 at finite length scales. Similarly to the 
bosonic case, 

z = l-^ + Oil/N^l (16) 

and the charge is marginally irrelevant to the order in 
1/N 0, M 

Besides the trivial fully attractive fixed point, the 
large-TV /3-functions in Eqs. (13)-(15) exhibit two crit- 
ical points in the attractive plane ^ gf = ga = cjx ~ 0: 
1) at ga = -1/4, 5c = and 2) gc = -1/4, ga = 0. 
There is also a bicritical point at ga ^ gc = ^1/4, which 
directs the flows towards one of the two critical points 
(Fig. 1). The critical points are related by the symmetry 
under a change of sign of 7^ for 'down' components with 
a = N/2 + 1, ...N accompanied by the exchange of ga and 
gc- The transition is either to A = cr^'cr^'cr) 7^ 0, 
which corresponds to an AF with a finite staggered mag- 
netization, or to a CDW, with the finite staggered density 
C = {J2a ^o-^o-) 0- The same, of course, follows from 
the explicit solution of the model at iV = 00. The flow of 
gj towards the origin also agrees with the saddle-point 
equations, which do not show a ferromagnetic critical 
point at N = 00, whereas the irrelevance of gd simply 
means that the chemical potential vanishes. 



Eqs. (14), however, appear to exhibit additional crit- 
ical points at g^. = 1/2. These, however, would occur 
within the AF or the CDW, and are artifacts of our pro- 
cedure which checks only the stability of the semi-metal. 
It is easy to see from the explicit solution that the ex- 
isting gap prevents such an additional transition. All g^ 
are therefore irrelevant. 

The transition in the pure, e = 0, repulsive Hubbard 
model with ga < and > in the large- iV limit is con- 
trolled therefore by the critical point A. Recalling that 
ga = Ua"^ A/ {8n), with Aw 1/a, and tay/i/2 = 1 by 
our convention, one finds that this corresponds to the 
critical value of U^/t ~ 5.5, certainly fortuitously close 
to the values found in numerical calculations 0. Above 
the critical interaction the system develops a gap, at the 
same time becoming insulating and antiferromagnetic. 
At the critical line gc = gf = Qd — gx = e^ — 0, upon 
the change of sign of 7-matrices for 'down' components 
the Lagrangian becomes identical to the much studied 
Gross-Neveu model in 2+1 dimensions [12. Evidently, 
the Gross-Neveu critical point has only one unstable di- 
rection to the order 1/A^. Since the actual expansion 
parameter is 4A'^, I expect this feature to survive even 
for N = 2. This leads to the conjecture that the SM - 
AF transition in the Hubbard model is continuous and 
described by the N = 2 Gross-Neveu critical point, at 
which 

{^M^)'^Aq,^)) - (g' + c.2)(''*-i)/2, (17) 

with the fermion's anomalous dimension 7711. = 
(2/(37r2/V)) + 0{1/N^) 17j. The order parameter's cor- 
relation function at the critical point also decays as: 

{A{x, t)A(0, 0)) ~ {x^ + r2)-(i+'')/2, (18) 

where 77 is the standard anomalous dimension, and 77 = 
1 - 16/{3tt'^N) + 0{1/N^). The correlation length di- 
verges at the critical point with the exponent v = 1 + 
8/(37r2/V) -t- 0{1/N'^), and the usual scaling laws are ex- 
pected to be satisfied. The critical exponents have been 
computed to the order 1 /N^ (with 77,1, known even to the 
order 1/N^) [l^, as well as being determined by Monte 
Carlo calculations, the e-expansion and the exact 
renormalization group [l9l |. In summary, for N — 2 
one finds 77^, = 0.038 ± 0.006, v = 0.97 ± 0.07, and 
77 = 0.770 ±0.016 [13. 

The presence of gapless fermions on the semi-metallic 
side places the Gross-Neveu phase transition outside the 
usual Ginzburg-Landau- Wilson paradigm, as evidenced 
by the large anomalous dimension 77, for example. In 
fact, the Gross-Neveu model probably defines the sim- 
plest such a universality class. Its distinct characteristic 
is the fermion's anomalous dimension 77^, which governs 
the disappearance of quasiparticles as the transition is 
approached on the semi-metallic side. Scaling dictates 
that the residue of the quasiparticle pole behaves as 

(t/c-t/)"*, (19) 
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so that a very small r|^J, would make it appear discontin- 
uous a.t U = Uc- 

In the AF, the eight generators that anticommute with 
cTz 7o become broken; these are {hjC^z) (X" (73, 75) and 
{ax,(7y) (g) (/, 735). Among these only {ax,cry) (K) /, which 
generate the usual spin rotations, correspond to the exact 
symmetry at J7 = 0, whereas the rest emerge as genera- 
tors of approximate symmetries only at low energies. In 
the insulating phase the Goldstone bosons which corre- 
spond to the emerging generators a re g apped, due to the 
irrelevant terms excluded from L |20j. The low-energy 
spectrum in the insulator consists therefore only of the 
usual magnons. 

The long-range nature of Coulomb interaction is found 
to be irrelevant at a large N. On the scale of lattice spac- 
ing, however. Coulomb interaction leaves its imprint on 
the initial value of the coupling gc, as indicated right be- 
low Eq. (11), for example 21]. In general, if the nearest- 
neighbor repulsion V is made sufficiently strong so that 
the line of initial conditions in Fig. 1 reaches left of the 
point C, there is an additional semi- metal - CDW tran- 
sition. Identifying V with ^ e^/a gives an alternative 
mechanism to that of ref. 8 for the CDW formation. 
The two lines of continuous transitions merge above a 



certain when the line of initial conditions comes left 
of the point B. The direct transition between the AF and 
the CDW is discontinuous. It seems natural, however, to 
assume that in reality U > V, which would suggest a 
single, continuous, antiferromagnetic transition. 

For graphene, t w 2.5eV, C/ w 5 - UeV, U/V w 2 - 3 
[lol | , so that the system is probably on the SM side of the 
transition. The external magnetic field, however, changes 
the density of states into a series of delta- functions, so 
that the transition can now in principle take place even 
at an infinitesimal coupling [Tll | . In the magnetic field the 
flow of the couplings should be cutoff at ^ , where 
3> a is the magnetic length. If the large-N picture 
presented here holds for N — 2, in the the pure, e = 0, 
Hubbard model, at a sufficiently low field all couplings 
would become negligible compared to ga- This would 
suggest that the magnetic field, at least with the Zeeman 
term neglected and the long-range component of the in- 
teraction screened by a metallic substrate, for example, 
should 'catalyze' the antiferromagnetic order at a weak 
U. 
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